We want to see if there are more stable loop orbits around certain axes using frequency mapping. I'll take a simplistic, 3 component Milky Way potential and run frequency mapping on 5 planes in x-z at the same energy but with different angles, $\phi$, in the x-y plane:



In [8]:

    
from __future__ import division, print_function

# Third-party
import astropy.units as u
import matplotlib.pyplot as plt
import numpy as np
%matplotlib inline

# Custom
import gary.potential as gp
from gary.units import galactic
from streammorphology.potential import potential_registry
from streammorphology.initialconditions import tube_grid_xz



In [15]:

    
E = -0.21
nplanes = 5
fig,axes = plt.subplots(1, nplanes, figsize=(nplanes*5,6), sharex=True, sharey=True)

for i,phi in enumerate(-np.linspace(0.,90.,nplanes)*u.deg):
    potential = gp.TriaxialMWPotential(phi=phi)
    w0 = tube_grid_xz(E, potential, dx=1., dz=1.)
    
    axes[i].plot(w0[:,0], w0[:,2], linestyle='none')

axes[0].set_xlim(-5,64)
axes[0].set_ylim(-5,64)
fig.tight_layout()



In [ ]: